Probing interaction in the dark sector 
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A phenomenological attempt at alleviating the so-called coincidence problem is to allow the dark 
matter and dark energy to interact. By assuming a coupled quintessence scenario characterized 
by an interaction parameter e, we investigate the precision in the measurements of the expansion 
rate H{z) required by future experiments in order to detect a possible deviation from the standard 
ACDM model (e = 0). We perform our analyses at two levels, namely: through Monte Carlo 
simulations based on eCDM models, in which H{z) samples with different accuracies are generated 
and through an analytic method that calculates the error propagation of e as a function of the error 
in H(z). We show that our analytical approach traces simulations accurately and find that to detect 
an interaction using H{z) data only, these must reach an accuracy better than 1%. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36. +x 



I. INTRODUCTION 

There is nowadays significant observational evidence 
that the Universe is currently undergoing accelerated ex- 
pansion. However, although fundamental to our under- 
standing of the Universe, several important issues about 
the mechanism behind cosmic acceleration as well as its 
role in the cosmic dynamics remain unanswered (see, e.g., 
for recent reviews). Among these open questions, the 
possibility of a non-minimal coupling between the two 
major energy components in the Universe, i.e. dark mat- 
ter and dark energy, has been widely investigated in the 
current literature [2]-[5]. 

Interacting dark matter/dark energy models violate 
adiabaticity and constitute a phenomenological attempt 
at alleviating the so-called coincidence problem (see, 
e.g., [6 ). In general, these models are characterized by a 
dilution of the dark matter density pm which is modified 
with respect to the usual scaling, i.e., 0^ 
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where e can take positive and negative values depending 
on if the transfer of energy is from dark energy to dark 
matter or vice versa, respectively. 

By introducing the above result into the balance equa- 
tion for the dark matter particles and dark energy, we 
find [5] 



From the observational point of view, current back- 
ground tests show that the interaction parameter e takes 
values close to zero [5 , which makes this class of models 
indistinguishable from usual ^CDM scenarios (e = 0). 
In principle, this degeneracy may be broken by studying 
the growth of structure in these models, which is modi- 
fied by the attenuated matter density evolution (Eq. [T]), 
and a possible distinction between coupling and uncou- 
pled models may be verified with the upcoming data from 
large redshift surveys [7]- 

In what concerns the background tests, it is worth 
mentioning that due to the multiple integrals that relate 
cosmological parameters to cosmological distances, direct 
determinations of the expansion rate place the tightest 
constraints on the dark energy equation of state w by 
reducing the so-called smearing effect [8 . In this paper, 
we investigate how measurements of the expansion rate 
H{z) may provide constraints on the interacting parame- 
ter e. To perform such analysis, we first use an analytical 
approach relating the accuracy of H{z) measurements to 
that of the interacting parameter e. Then, we generate 
H{z) samples from Monte Carlo simulations with differ- 
ent accuracy and values of the interaction scale e based 
on the matter evolution law ([T]) . By considering values of 
e that are consistent with current observations (e < 0.1), 
we show that a possible distinction between interacting 
and non-interacting models may be achieved only with 
an accuracy in H{z) data better than 1%. 
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II. ANALYTICAL CONSTRAINTS FROM H(Z) 



where stands for the dark energy density and we have 
assumed the equation-of-state parameter w = P(f)/P(f) re- 
lating p0 and the dark energy pressure p^. 
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We restrict the present analysis to coupled 
quintessence models with w = —1, which is math- 
ematically equivalent to the so-called dynamical A 
scenarios. In this case, Eq. Q reduces to 
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where the subscript denotes present-day quantities. As 
discussed in Ref. [4, such a evolution law seems to be 
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FIG. 1: Analytic computation of the relative error in H(z) as a function of the relative error in e for z = 0.5, 1.0 and 1.5. The 
left and middle panels show the case e = 0.1 and e — 0.05, respectively. The redshift dependence of ^ for values of e = 0.1 
and ^ = 20%, 50% and 80% is shown in the right panel. 



the most general one, having many of the previous phe- 
nomenological attempts as a particular case (see, e.g.. 
Table I of [9 ). 

The Friedmann equation for this kind of models is 
given by 



8H over a redshift interval by using 
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where the baryon, CDM and dark energy densities are 
given in units of the present critical density, x = 1 -\- z 
and Q\ci = - As written above, Eq. fl ex- 
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presses the fact that we are considering a non-mmimal 
coupling between the dark energy component and CDM 
particles only. We refer the reader to Ref. [10 for a dis- 
cussion about bounds on a possible interaction between 
dark energy and conventional matter from local gravity 
experiments and from primordial nucleosynthesis. 

From the above equation, we calculate the relative er- 
ror in the expansion rate as a function of the relative er- 
ror in the interaction parameter e from SH^ = (^) Se'^. 
This is shown in Fig. [ijfor values of e = 0.1 (left panel) 
for z = 0.5, 1.0 and 1.5. Note that a 0.5 (1)% error in 
H dit z = 1.0 results in a 40 (80)% uncertainty in e, but 
the precision varies with redshift (right panel). For the 
sake of comparison, we show in the middle panel of Fig. 
[l]the analysis for e = 0.05. As expected, the smaller the 
value of the interacting parameter the better the accu- 
racy in H{z) measurements required to distinguish be- 
tween this kind of scenario and non-interacting models 
with e = 0. In the right panel of Fig. 1 we show the 
dependence of 5H/H with redshift for the case e = 0.1. 
For all the analysis performed in this paper, we consider 
a flat universe with Hq = 73 km.s~^.Mpc~^, Q^o = 0.04 
and ftmo = 0.23 fTTl. 

The above analytical approach poses difficulty when 
one tries to compare these predictions with data. The 
reason is that it depends on a specific value of the red- 
shift whereas to fit the data (real or simulated) we use 
information of various data points spread over a redshift 
interval. In order to take that into account we averaged 
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where Zi = 0.1 and Zf = 1.5. 

The result is shown in Fig. [5] (left) for three different 
values of e. Note that the relative errors in e increases 
considerably with Therefore, we restrict ourselves to 
analyze the accuracy required in H{z) measurements to 
distinguish between a eCDM model and the standard cos- 
mological scenario (e = 0), which is formally determined 
at the point 5e = e (100% error). According to these re- 
sults, we see that if e = 0.1, an accuracy of SH/H < 1% 
would be necessary to rule out ACDM with la confidence 
level. For a weaker interaction, however, e.g. e = 0.05, 
SH/H must be smaller than 0.5%. 



III. SIMULATIONS OF H(Z) 

Current data of the Hubble parameter has roughly 15% 
of uncertainties [12 . At such a level, it is not possible 
to distinguish between the standard ACDM model and 
a interacting scenario with e < 0.1. On the other hand, 
future surveys will be able to measure the expansion rate 
H{z) with better precision. This is the case, for instance, 
of the Southern African Large Telescope [13] and the At- 
acama Cosmology Telescope [14 , which may reach 3%- 
10% accuracy from age estimates of passively evolving 
galaxies at high-z and large photometric surveys, like J- 
PAS, that may reach a 2% level from measurements of 
the radial scale of the baryonic acoustic oscllations [7] . 

To anticipate the aforementioned project's results con- 
cerning a possible interaction in the dark sector, we made 
use of Monte Carlo simulations to generate samples of 
H{z) data. The fiducial model in our simulations is a 
eCDM cosmology, represented by Eq. Q, with the set 
of parameters P = (i^o, ^60? ^mo) mentioned earlier and 
some selected values of the interaction parameter e in the 
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FIG. 2: Right: Similar to Fig. [T] but averaging over the redshift interval [0.1-1.5]. Left: A Monte Carlo realization of 15 data 
points of the Hubble parameter with 1% accuracy. 



ranse [-0.2 -0.2]^ For each value of e, we simulated sam- 
ples with relative error in H{z) of 3%, 1%, 0.5%, 0.3% 
and 0.1%. For each of these combinations, we follow 
Ref. [T5 and run 500 simulations of 15 points equally 
spaced in redshift corresponding to observations in the 
range [0.1 - 1.5]. 

The simulated data of H{z) are drawn from a normal 
distribution N^ji^cr) with mean ji = Hfid (see Eq. [4| 
and standard deviation a = (^) i^/id, where (^) is 
the relative error in the determination of H{z). In or- 
der to calculate the error bars we made a bootstrapping 
study of the current observational errors of H{z) data, as 
obtained by Stern et al. p!2] using age difference between 
passively evolving galaxies at different and estimated 
the deviation S from the average 15% value. Suppos- 
ing the errors have a normal distribution, the error bars 
are drawn from the normal distribution gh = N{cf^ Scr). 
Fig. 2 (left) shows a Monte Carlo realization of 15 simu- 
lated values of the Hubble parameter with 1% accuracy 
by assuming e = 0.1. 



IV. RESULTS 

Consider a fiducial model characterized by a given 
value of e = and the set of parameters P given above. 
For each of the 500 realizations, we minimize the function 
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The average e and its standard deviation are determined 
from the distribution of the resulting 500 best-fit values of 
e. Figure [3] shows the main results of our analyses. Plots 



^ Although negative values of e are forbidden by thermodynam- 
ics 21 3 observational data do not exclude this possibility |5j. 



of e with error bars corresponding to la standard devi- 
ation from the mean are shown for the cases of e = 0.1 
(left) and e = 0.05 (right). From the left panel, we clearly 
see that if the Hubble parameter reaches a relative error 
of ^ < 1%, we would be able to distinguish between a 
eCDM model with e = 0.1 and the standard ACDM sce- 
nario at Icr level. For e = 0.05 (right panel), a distintion 
at 1(7 level between these two classes of models is possible 
only when the relative error of H{z) observation reaches 
^ < 0.5%. 

For the sake of comparison, we also plot in Fig. |3]the 
analytical result derived from Eq. ([5|. Note that it (blue 
line) shows almost no difference from the simulations, 
proving to be an accurate tool for cosmological forecasts. 
It is worth mentioning that we also performed the same 
analyses presented in Fig. [3] by considering negative val- 
ues of e. In this case, we note that this range of values 
requires a slightly less accuracy in H{z) measurements, 
with ^ ranging from ~ 5 — 10% higher than those corre- 
sponding to positive values. The theoretical results from 
Eq. ([5| and simulations also agree for values of e < 0. 



V. CONCLUSIONS 

We performed Monte Carlo simulations of the expan- 
sion rate H{z) in a universe with interaction in the dark 
sector. The analyses have been performed by consider- 
ing values of the interacting parameter e in the interval 
[-0.2-0.2], which encompasses the current observational 
bounds on this quantity and have assumed uncertainties 
on H{z) of ^ = 3%, 1%, 0.5%, 0.3% and 0.1%. For each 
combination of ^ and e we generated 500 Monte Carlo 
simulations of H{z) in the redshift range [0.1 - 1.5]. 

We used these simulations to investigate the possibility 
of future projects probing the interaction between dark 
matter and dark energy through the model described by 
Eq. Q. We found that values of e that are within cur- 
rent estimated errors (0.1 and 0.05) might be detected 
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FIG. 3: The average fitted value of e as a function of tiie relative error in H(z). The left and right panels show, respectively, 
the fiducial values (dashed line) e = 0.1 and e = 0.05. The blue curves stand for the analytical result derived from Eq. 



through H{z) observations only if SH/H < 1%. It should 
be noted that the observational accuracy can be appre- 
ciably increased by combining H{z) measurements with 
data from other observables, such as SNe la and CMB 
data. Currently, these observations weakly constrain e 
and in some cases favor negative values [5 . 

We also discussed an analytical procedure to foresee 
5e/e as a function of SH/H [Eq. ([s])] and showed that it 
provides results in agreement with Monte Carlo simula- 
tions. We emphasize that this procedure can be applied 
to any observable as well as to joint analyses involving 
different cosmological probes. We plan to extend this 



analytical approach to combined probes in a forthcom- 
ing communication. 
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